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STOCHASTICITY  OF  PARTICLE  ORBITS 
IN  A  COLLECTIVE  PARTICLE  ACCELERATOR 

I.  INTRODUCTION 


High  energy  and  high  density  beams  have  been  of  considerable  interest 
in  recent  years.  A  novel  method  to  achieve  such  a  high  density  high  energy 
particle  beam  was  proposed  by  Friedman1  in  1978  in  a  device  called  the 
Collective  Particle  Accelerator  (CPA).  Briefly,  the  device  is  a  travelling 
wave  accelerator  with  a  novel  feature:  the  travelling  wave  (T.W.)  which 
does  the  accelerating  is  generated  on  the  axis  of  a  cylindrical  drift  tube 
as  a  consequence  of  the  interaction  between  (1)  a  periodically  bunched 
annular  electron  beam  injected  near  the  drift  tube  wall,  and  (2)  a 
periodically  rippled  axial  magnetic  field.  The  advantage  of  this  as 
compared  to  conventional  travelling  wave  devices,  in  which  the  T.W.  is  a 
cavity  mode,  arises  from  the  fact  that  the  phase  velocity  in  the  CPA  may  be 
chosen  to  be  as  slow  or  fast  as  desired,  within  limits  to  be  discussed 
below,  and  may  be  easily  varied  as  a  function  of  distance  along  the  drift 
tube  by  varying  the  wavelength  of  the  magnetic  ripple.  Thus  ions  may,  in 
principle,  be  accelerated  along  the  axis  either  coparallel  or  antiparallel 
to  the  motion  of  the  annular  beam,  or  electrons  may  be  accelerated  in  the 
antiparallel  direction.  An  important  side  effect  of  this  method  of 
generating  the  T.W.  is  the  fact  that  three  large  amplitude  waves  are 
simultaneously  generated  in  the  device:  a  forward  (i.e.  coparallel  with  the 
annular  beam  velocity)  T.W.;  a  backward  T.W.;  and  a  stationary  wave 
pattern. 

The  first  experiments  in  support  of  the  CPA  concept  were  conducted 
around  1976.  A  relativistic  annular  electron  beam  was  propagated  a 

distance  of  7  meters,  without  bunching  and  with  straight  magnetic  field 
lines.  In  the  next  phase  the  beam  was  bunched  and  again  7  meter 
propagation  was  achieved.  Addition  of  the  magnetic  field  ripple  in  the 
third  stage  of  the  effort  was  similarly  successful.  Next,  a  solid  beam 
with  current  -10%  of  the  annular  beam  current  was  successfully  propagated 
along  the  axis  in  the  absence  of  the  annular  beam.  The  full  experiment, 
solid  beam  plus  bunched  annular  beam  with  rippled  B  field  was  then  tried. 
The  solid  beam  did  not  come  out  at  the  end  opposite  to  which  it  was 
injected,  and  diagnosis  of  this  beam  through  the  annular  beam  was 
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difficult,  so  a  simpler  experiment  was  tried.  A  solid  beam  was  injected 
into  an  annular  beam  without  bunched  or  B  field  ripples,  but  was  not 
detected  near  the  other  end,  although  a  magnetic  field  disturbance 
attributed  to  the  solid  beam  was  detected  close  to  the  solid  beam 
emitter.  Note  that  without  bunches  and  ripples  no  acceleration  can  take 
place,  in  which  case  the  solid  beam  cannot  overcome  the  potential  hill  near 
the  annular  beam  cathode,  so  in  the  full  experiment  one  does  not  expect  the 
solid  beam  to  hit  the  opposite  end  even  if  it  propagated.  A  number  of 
variations  of  these  experiments  were  attempted  in  1979;  e.g.,  the  solid 
beam  was  carried  to  a  point  well  inside  the  drift  tube  via  a  smaller  radius 
drift  tube  projecting  inside  the  annular  beam,  but  still  the  solid  beam  was 
not  observed  to  propagate  the  remaining  distance.  Again  we  emphasize  that 
observation  of  the  weak  solid  beam  through  the  strong  annular  beam  is 
difficult,  and  the  negative  result  in  the  unmodulated  case  cannot  be 
regarded  as  definitive. 

Various  theoretical  analyses  for  the  stability  properties  of  the  dual 
beam  system  were  undertaken  in  order  to  understand  the  experimental 
results.  Chen  and  Palmadesso2  studied  the  diocotron  instability  in  the  CPA 
in  an  attempt  to  understand  the  sensitivity  to  small  misalignments  of  the 
annular  beam  encountered  in  the  first  experiment.  The  l  *  1  diocotron 
mode  can  be  excited  by  a  shift  of  the  annular  beam  axi3  off  the  drift  tube 
axis,  and  was  shown  to  have  a  growth  length  cmparable  to  the  device  length 
for  the  parameters  used.  The  growth  of  the  instability  is  not  significant 
if  the  beam  is  very  accurately  centered  initally,  but  is  enough  to  cause 
the  beam  to  hit  the  wall  if  the  initial  displacement  is  within  an  e-fold  or 
two  of  the  distance  between  the  beam  and  the  wall. 

Uhm^  showed  that  the  transverse  oscillation  instability  has  a  growth 
time  -  1/5  of  the  beam  transit  time  for  the  CPA  parameters  used.  However, 
Ganguli  and  Palmadesso4  used  the  wave  kinetic  equations  in  a  finite  length 
device  to  show  that  the  transverse  oscillation  instability  behaves  like  a 
backward  convective  instability  in  the  rest  frame  of  the  solid  beam  i.e., 
in  this  frame  the  growth  is  spatial  rather  than  temporal  and  vanishes  at 
the  beam  head.  The  implication  is  that  the  front  part  of  the  beam  (about  5 
meters  worth  in  this  case)  would  make  it  across  the  device  even  though  the 
tail  end  of  the  beam  was  being  torn  up  as  the  head  completed  its  transit. 
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Thus  the  theoretical  models  analyzed  so  far  do  not  account  for  the 
failure  of  the  solid  beam  to  emerge  from  the  end  of  device,  in  the  case 
where  an  accelerating  T.W.  is  present.  These  theoretical  models  did  not 
investigate  the  role  of  the  important  second  and  third  large  amplitude 

waves  generated  and  the  influence  of  these  waves  on  the  particle 
dynamics.  The  purpose  of  this  paper  is  to  study  the  role  of  these  waves  in 
determining  the  particle  dynamics  of  the  CPA  and  hence  influencing  the 

experimental  outcome.  In  section  II  we  describe  our  theoretical  model, 
provide  the  equation  of  motion  governing  the  CPA  system  and  describe  the 
method  of  solution.  In  section  III  we  give  the  results  of  the  calculations 
and  discuss  the  ramifications  of  these  and  finally  in  section  IV  we  provide 
concluding  remarks. 

II.  THEORY 

As  discussed  earlier  an  important  side  effect  of  the  way  that  the 
travelling  waves  are  generated  in  the  CPA  is  the  simultaneous  presence  of 
three  large  amplitude  waves  in  the  device.  It  is  well  known  that  if  a 
particle  is  under  the  influence  of  two  or  more  waves  with  resonance  widths 
broad  enough  to  bridge  (or  nearly  bridge)  the  gap  between  the  wave  phase 

velocities,  then  the  multiple  resonance  can  make  the  particle  orbits 

stochastic5.  With  the  onset  of  stochastic ity  the  predictability  of  the 
particle  motion  is  lost.  This  is  a  possible  scenario  taking  place  inside 
the  CPA.  The  expression  for  the  axial  electric  fields  generated  by  the 
bunched  beam  and  rippled  B  field  131 , 
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where  L  is  the  wavelength  of  the  rippled  magnetic  field,  f  and  X  are  the 

frequency  and  the  wavelength  of  the  modulation.  E  -  (3/4e  )  (r,/r  L) 

zo  o  1  0 

where  r  ,  r1  are  the  equilibrium  radius  of  the  Annular  Electron  Beam  (AEB) 
and  the  amplitude  of  the  oscillations  of  the  AEB  due  to  the  ripples  of  the 
magnetic  field,  and  Q  is  charge/length.  Equation  (1)  represents  the 
superposition  of  a  forward  wave  with  a  phase  velocity  given  by, 
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where  v  -  fx,  a  backward  wave  with  a  phase  velocity  of 


and  a  stationary  wave  pattern. 

Both  travelling  waves  can  accelerate  ions  but  only  the  backward  wave  can 
accelerate  the  electrons  since  it  can  have  a  phase  velocity  approaching  c 
(the  velocity  of  light)  which  happens  to  be  approximately  the  velocity  of 
the  electrons  parallel  to  the  magnetic  field. 

A  particle  subjected  to  the  field  given  in  (1  )  experiences  the 
potential  of  both  the  forward  and  the  backward  waves  along  with  the 
standing  wave  simultaneously.  These  are  large  amplitude  waves  whose 
resonance  widths  in  velocity  space  scale  as, 

ia  .  i  (fsj"2  .  (rm  >-’'2  <4) 
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o 

where  Aoj  is  the  resonance  width  in  the  frequency  space,  k  is  the  wave 

2  2  -1/2 

number  and  Y  »  (1  -  V  /c  )  .  For  the  typical  CPA  parameters 

1  /2 

(eE  k/m  )  /k  -  0(c)  ;  for  low  Y  this  implies  an  overlap  of  resonance 
zo  o 

leading  to  stochastic  motion.  As  a  consequence  of  the  stochasticity  in  the 
orbits,  the  particles  (electrons  in  our  case)  spend  longer  time  inside  the 
device  thereby  increasing  the  space  charge  and  decreasing  the  value  of  the 
critical  current.  However,  equation  (^)  also  implies  that  for  sufficiently 
large  injection  energy  (initial  value  of  Y)  the  resonance  widths  will  be 
small  enough  so  that  stochastic  motions  will  not  occur.  In  view  of  the 
results  of  the  CPA  experiments  we  must  investigate  the  role  of 
3tochasticity  and  estimate  injection  energy  thresholds. 


The  equation  of  motion  for  the  electrons  in  the  CPA  is  given  by, 
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where  e  is  the  magnitude  of  electronic  charge  and  E_  given  by  (1 )  can  be 
expressed  as, 
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Here,  c  -  z/X,  a  -  X/L  and  t  •  ft  +  1/4  .  The  equations  of  motion  can  be 
rewritten  in  dimensionless  form  as, 
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where  a  »  c/fX,  p  -  p/fXmo>  mQ  is  the  rest  mass  of  the  electrons. 

We  shall  first  study  an  infinitly  long  system  by  means  of  surface-of- 
section  plots.  After  injecting  a  single  particle  at  t  =  0  we  shall  let  the 
system  represented  by  (8)  and  (9)  run  in  time,  and  plot  p  and  5  for  x  => 
1,2,3....  etc.,  given  the  initial  momentum  p  and  position  t  .  We  take  E_0 
»  15  MV/meters,  a  =  2,  f  »  .995  x  10y  Hz,  and  X  -  0.3  meter.  These 

parameters  are  typical  of  the  experiment1.  Note  that  the  system  is 
periodic  in  q  and  hence  we  shall  plot  5  modulo  1  .  The  results  will  be 


discussed  in  the  following  section 


After  examining  the  infinite  device  to  gain  general  insights  we  will 
study  a  finite  device  of  7  meters.  Accumulation  of  electrons  in  the  device 
is  detrimental  to  its  operation  and  may  be  caused  by  delays  in  time  of 
flight.  This  enhances  the  electron  density  inside  the  device  and  decreases 
the  critical  current.  In  order  to  estimate  the  lag  time  of  the  electrons 
we  plot, 


as  a  function  of  time.  Here  t.  is  the  transit  time  of  the  electrons 
through  the  device  and  =*  7/c  is  the  minimum  time  that  the  electron 

needs  to  travel  the  device  of  7  meters  when  the  electron  velocity  remains 
close  to  the  speed  of  light  c.  Thus  orbits  with  small  S  do  not  contribute 
to  charge  accumulation  in  the  device  while  the  higher  values  of  S  indicate 
highly  distrubed  electron  orbits  which  contribute  to  an  enhancement  of  the 
space  charge  cloud  and  a  decrease  in  the  limiting  current. 

Also,  experimentally  it  is  not  possible  to  launch  all  the  particles  at 
the  same  phase  initially.  We  shall  study  S  as  a  function  of 
$  »  -2ir(^o  +  tq  -  1/4),  the  value  at  injection  of  the  backward 

(accelerating)  wave  phase,  for  various  injection  energies  Y  .  The  results 
are  discussed  in  the  following  section. 

III.  RESULTS  AND  DISCUSSION 

We  find  two  kinds  of  orbits  exhibited  on  the  surface  of  section 
plot.  "Integrable"  orbits  are  indicative  of  the  existence  of  a  well- 
behaved  constant  of  the  motion  which  confines  the  motion  of  the  particle 
through  phase  space  to  a  subspace  of  reduced  dimension.  The  intersection 
of  an  integrable  orbit  with  a  two-dimensional  surface  of  section  such  as 
that  in  Fig.  1  takes  the  form  of  a  sequence  of  points  [p(xn),c(T  )] 
for  in  »  0,1,2,...,  which  lie  on  a  simple  one-dimensional  curve  or,  in 
special  cases,  on  a  single  "fixed"  point.  In  general,  the  obvious  global 
constants  of  the  motion  are  exploited  to  reduce  the  dimension  of  the  system 
of  motion  equations  before  one  constructs  a  surface  of  section  plot,  so 
that  the  additional  constants  associated  with  integrable  orbits  exist  only 


in  local  regions  of  phase  space,  and  it  is  rarely  possible  to  express  them 
analytically.  Each  of  the  nested  closed  curves  near 

5  mod  1  -  .5,  Y  <  40  in  figures  1  and  2  represents  an  integrable  orbit 

corresponding  to  a  different  initial  condition.  At  the  center  of  this 

group  of  nested  orbits  is  a  fixed  point,  which  corresponds  to  periodic 
motion  passing  through  the  stable  force  null  of  the  backward  wave  at 
integer  values  of  t  ,  i.e.,  the  electron  returns  to  the  point 

5  »  .5  mod  1  at  i  *  1,2, 3, etc.. 

The  four  or  five  outermost  curves  which  pass  through  5  mod  1  =  0 
with  p  <  -10  (i.e.  Y  >  10)  also  represent  distinct  integrable  orbits.  No 
calculations  were  made  in  the  blank  regions  at  higher  values  of  Y ,  but  it 
is  believed  that  all  orbits  in  this  region  are  integrable.  At  high 
energies  the  resonance  widths  of  the  three  waves  are  small,  the  resonances 
are  well  separated,  and  3olid  beam  electrons  (v  -  -c)  interact  strongly 

only  with  the  backward  wave.  The  electron  velocity  exceeds  the  backward 
wave  phase  speed  when  Y  >  10,  so  that  an  accelerated  electron  cannot 

maintain  a  constant  phase  relationship  with  the  wave.  Electrons  gain 

energy,  slip  smoothly  out  of  the  phase  interval  in  which  the  force  due  to 
the  electric  field  is  coparallel  with  the  velocity,  then  decelerate  while 
continuing  to  advance  in  phase.  As  the  particle  loses  energy  the 

resonances  broaden.  If  the  electron  in  this  idealized  infinite  system  has 
sufficient  energy,  it  will  slide  forward  in  phase  to  the  next  accelerating 
region  before  its  speed  decreases  to  the  point  where  it  matches  the  phase 
speed  of  the  backward  T.W.,  and  before  resonance  overlap  occurs.  The 
lowest  energy  integrable  orbit  of  this  group  of  orbits  slows  down  until  its 
speed  is  infinitesimally  greater  than  the  backward  wave  phase  velocity  at 
the  point  it  reaches  the  unstable  (hyperbolic)  force  null  of  the  backward 
T.W.  at  2ir(c+x)  -  0  (5  mod  1  -  0  on  the  surface  of  section  plots);  thus  Y 
never  falls  below  10  on  this  orbit.  There  is  also  a  set  of  integrable 

orbits  representing  high  energy  particles  with  positive  velocities  which 
are  outside  the  resonaces  of  all  three  waves.  These  orbits  are  not 

physically  meaningful  for  the  CPA  and  will  not  be  discussed  further  here. 

"Stochastic"  orbits  occupy  regions  in  which  no  such  local  well  behaved 

constants  of  the  motion  exist,  and  these  orbits  meander  through  the  entire 

phase  space  allowed  by  the  global  constants  of  the  motion,  with  the 

restriction  that  they  may  not  cross  the  "KAM"  surfaces^  (roughly  speaking, 


these  are  surfaces  in  phase  space  traced  out  by  integrable  orbits).  In  the 
2-D  surface  of  section  plots,  a  single  stochastic  orbit  accounts  for  all 
the  points  between  the  integrable  regions  discussed  above.  This  orbit 
behaves  in  a  chaotic  manner  at  low  energies  (Y  less  than  about  10),  but 
catches  the  accelerating  phase  of  the  backward  wave  at  random  intervals  and 
then  gains  energy  in  an  orderly  way.  Segments  of  the  stochastic  orbit  with 

Y  above  about  10  rise  and  fall  in  energy  in  much  the  same  manner  as  the 
high-energy  integrable  orbits,  but  the  electron  speed  falls  below  the 
backward  T.W.  phase  speed  (Y  -  10)  while  the  particle  is  still  in  the 
deceleration  phase,  so  that  the  particle  continues  to  lose  speed  and  energy 
and  slips  in  phase  with  respect  to  the  backward  T.W.  at  an  increasing 
rate.  Finally  the  motion  ceases  to  be  dominated  by  the  backward  T.W.,  and 
the  particle  falls  back  into  the  chaotic  resonance  overlap  region  below 

Y  -  10. 

A  rough  calculation  of  the  electron  orbits  for  Y  >  10  may  be  performed 
by  assuming  that  the  particles  have  a  constant  velocity  of  -c  and  are  only 
affected  by  the  backward  wave,  in  which  case  the  momentum  equation  takes 
the  form 

eE 

_£  ,  [ — |2_J  C03  [  >2ir(co+  (1-o)t  -  1/4)  ],  (11) 

if  m 

o 

where  is  the  injection  value  of  q.  This  equation  may  be  solved  to  give 
an  analytic  expression  for  the  energy  Y  of  the  particle  as  a  function  of 
the  phase  <p  of  the  particle  relative  to  the  backward  wave,  the  injection 
energy  y  ,  and  the  injection  phase  $: 

Y  '  Yo  +  fei2)  2^T)  [  sin(*}  “  sin(*}  )  (12) 

o 

-  Yq  +  280  (  sin($)  -  sin(ij))  ) 

where  the  CPA  parameters  have  been  used.  An  electron  injected  at  the 
unstable  force  null  ($  =  n/2)  at  the  speed  of  the  backward  wave  (yq=  10' 
achieves  its  maximum  energy  at  phase  <p  =  — ir/2.  This  energy  is  Y  =  570, 
which  is  in  good  agreement  with  Fig.  1  for  the  threshold  energy  above  which 
all  electrons  have  integrable  orbits.  In  a  finite  length  device,  the 
electron  stays  with  the  backward  wave  only  for  a  finite  time  and  the  phase 


slippage  is  limited  in  size.  Thus,  not  all  particles  below  the  threshold 

will  become  stochastic  in  a  finite  length  device.  For  a  device  of  length  7 

meters,  a  particle  traveling  at  the  speed  of  light  will  remain  in  the 
-8 

device  for  2.3x10  sec,  which  implies  that  At  -  23  for  this  choice  of  the 
parameters  f  and  X,  and  that  the  maximum  phase  shift  that  the  particle  can 
undergo  is  2it(1-<j)At  -  -ir/4  .  From  equation  (12),  this  means  that, 
if  Y  >  1 0  +  280(sin(ir/8)  -  sin(-ir/8))  -  224,  a  particle  cannot  become 
stochastic  in  a  7  meter  device  with  these  parameters. 

In  figures  3a  to  3g  we  display  plots  of  the  normalized  lag  time  S 
against  the  injection  phase  of  4>  of  the  backward  wave  for  7  different 
injection  energies  from  Y  -  5  to  Y  =■  250.  These  plots  are  pertinent  to  a 
finite  device  of  7  meters.  We  learn  two  very  important  features  of  the  CPA 
operation  from  this  figure.  First  and  foremost  it  is  seen  that  as  the 
injection  energy  Y  increases  from  Y  -  5  to  Y  *  250  the  mean  value  of  S 
decreases  steadily.  This  implies  that  for  high  enough  injection  energies 
there  is  uniform  motion  no  matter  what  the  initial  phase  is.  For  lower 
energies  the  mean  value  of  S  increases  significantly  indicating  that  the 
electrons  are  spending  longer  times  inside  the  device.  This  implies  an 
electron  space  charge  cloud  build  up  inside  the  device  and  consequently  a 
decrease  in  the  limiting  current.  The  value  of  the  injection  energy  in  the 
experiment  for  the  solid  beam  was  on  the  order  of  Y  -  5,  and  hence  it  is 
not  surprising  to  find  the  absence  of  the  accelerated  solid  beam  in  the 
experiment,  which  was  conducted  with  injected  current  very  close  to  the 
critical  current  as  traditionally  calculated. 

Along  with  the  decrease  in  the  mean  value  of  S,  we  also  notice  a 
shrinkage  of  the  interval  A$  over  which  S  is  significant  as  the  injection 
energy  Yq  is  increased.  If  this  interval  of  injection  phases  is  avoided, 
then  uniform  motion  could  be  achieved  even  at  lower  energies.  For  example, 
if  Yq  *  5  and  $  =*  it,  figure  3a  shows  no  significant  lag  time.  On  the  other 
hand  even  for  high  energy,  3ay  Y  »  150,  (figure  3e)  if  the  initial  phase  is 
around  zero  we  can  expect  stochastic  behavior.  As  noted  earlier  we  find 
that  the  optimum  phase  for  acceleration  is  $  -  tt.  Figure  3a  indicates  that 
particles  launched  with  an  injection  phase  $  ■  n  can  be  accelerated  even 
for  injection  energies  lower  than  the  Y  of  the  background  wave  since  the 
particles  cannot  shift  into  the  chaotic  region  before  leaving  the  finite 


device.  When  a  uniformly  solid  beam  is  injected,  particles  of  all  phases 
will  simultaneously  occupy  the  drift  tube.  Thus  the  charge  build  up 

associated  with  the  particles  with  phase  in  the  stochastic  range  can  be 
detrimental  to  successful  acceleration  of  those  particles  which  have  the 
optimum  phase. 

IV.  CONCLUSION 

We  have  studied  the  influence  of  the  three  large  amplitude  waves  that 
arise  due  to  the  CPA  wave  generation  scheme  on  the  particle  dynamics  in  the 
CPA.  For  typical  CPA  parameters  the  resonance  width  of  the  waves  are  broad 
enough  to  overlap  for  low  energy  electrons.  This  gives  rise  to 

stochasticity  in  the  particle  orbits  if  the  initial  phase  is  not  properly 
chosen.  The  stochasticity  enhances  the  transit  time  of  the  electrons 
consequently  giving  rise  to  an  electron  cloud  build  up  inside  the  device. 
The  charge  build  up  inside  the  device  decreases  the  critical  current  and 
could  prevent  the  propagation  of  the  solid  beam.  Further,  we  conclude  that 
in  the  CPA  experiment  those  particles  of  the  solid  beam  that  were  launched 
at  the  optimum  phase  must  not  have  caused  any  problems.  Since  the  beam  was 
solid  there  was  a  continuous  range  of  injection  phases  and  the  particles  in 
the  stochastic  range  contributed  to  the  charge  build-up  in  the  device  which 
in  turn  prevented  propagation  and  acceleration  of  the  optimum  phase  locked 
particles.  To  avoid  these  effects  for  a  solid  beam  an  injection  Y  i  225 
is  necessary.  However,  if  the  solid  beam  is  injected  in  small  pulses 
around  the  optimum  phase  it  is  possible  to  achieve  successful  acceleration 
even  for  low  injection  energies. 

We  have  studied  the  lag  time  as  a  function  of  the  phase  of  the 
backward  wave.  It  seems  that  by  properly  choosing  the  injection  phases 
(pulsed  solid  beam)  it  i3  possible  to  avoid  the  stochasticity  and  minimize 
the  transit  time. 
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Fig.  3.  A  plot  of  the  lag  time  S  against  the  injection  phase  $  of  the 
backward  wave  for  a  device  of  7  meters.  The  initial  values  of  Y 


at  injection  are: 

(a)  Y  -  5,  (b)  Y  -  10,  (c)  Y  -  15,  (d)  Y  -  50,  (e)  Y  -  150,  j 

(f)  Y  -  200  ,  and  (g)  Y  -  250.  - 
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